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Observations made during the MIZEX program indicate the presence of 
mesoscale eddies in the ocean front at the marginal ice edge in the East Greenland 
Current. The eddies ranged in scale from 5 to 80 km. Barotropic and baroclinic 
instability may be the physical mechanisms responsible for the existance of such eddies. 
The observations also indicate transient wind reversals (3-10 m/s) with a frequency of 
several days. Here the effect of time-dependent winds and ocean eddies on ice motion 
in a marginal ice zone is studied. Results are obtained w r ith a two-layer, nonlinear, 
primitive equation ocean model and a coupled free-drift ice model. The results indicate 
that ocean eddy signature in the ice edge is sensitive to cross-ice-edge motion induced 
by the winds and is shown to be dependent on magnitude, direction, and duration of 
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I. INTRODUCTION 



A. THE MARGINAL ICE ZONE AND THE EAST GREENLAND CURRENT 

The marginal ice zone is a region where complex atmospheric and oceanic 
physical processes afreet the distribution of ice between the pack ice and the open 

ocean. In recent years the marginal ice zone (MIZ) in the region of the East 

Greenland Current has been extensively investigated during Marginal Ice Zone 

Experiments. The East Greenland Current is arctic (polar) water flowing southward 
along the east coast of Greenland, over the 200 m deep and 200 km wide continental 
shelf, and then continuing to flow southward over deeper (4000 m) water (Wadhams et 
al . , 1979). The boundary formed by the arctic water of the East Greenland Current 
and the Atlantic water is a potentially unstable ocean front. Figure 1.1 shows the 
location of the East Greenland Current and the mean ice edge. 

The study of stability of ocean fronts has been done in the laboratory as well as 
in numerical experiments. Griffiths and Linden (1981) examined stability 

characteristics of density driven fronts in rotating tank experiments. Their experiments 
indicate that the type of instability (barotropic or baroclinic) associated with a given 
flow field is dependent upon the length scale. In general, if the length scale L of the 
mean flow is larger than the local internal Rossby radius of deformation (Rj), then the 
resulting instability is baroclinic. Where the length scale is comparable to R^, the 
instability is more likely to be barotropic. The width scale of the East Greenland 
Current is on the order of 200 km, much greater than the local Rj (8 km), suggesting 
baroclinic instablitv (Griffiths and Linden, 1981). The observed 50 km length scales of 
some of the eddies (Wadhams and Squire, 1983) appear to be consistent with these 
considerations. The mesoscale ice edge frontal eddies observed by Johannessen et al. 
(1983) had a length scale comparable to the local internal Rossby radius; therefore, 
from the above considerations, these mesoscale eddies would be the result of barotropic 
instability. 

B. EDDIES IN THE EAST GREENLAND CURRENT REGION 

From an observational standpoint, eddies in the MIZ have only recently been 
documented. These eddies have been observed to have scales ranging from 5 to 80 km. 
There are five possible sources for the eddies of the East Greenland Current region: 
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Figure 1.1 Ice-Edge Oceanic Front (from Johannessen et al., 1987). 

(1) Barotropic instability was suggested by Johannessen et al. (1983) as being 
responsible for the formation of small-scale eddies; (2) Baroclinic instability is largely 
responsible for the larger-scale eddies (Griffiths and Linden, 1981); (3) Current- 
topography interaction may cause eddy formation (Smith et al ., 1984); (4) Open ocean 
eddies present in Atlantic water may be advected toward the ice edge, leading to 
interaction that can develop into ice-edge eddies (Johannessen et al., 1987); and (5) Ice 
edge air-sea interaction may cause eddies to form (Roed and O'Brien, 1983; 
Hakkinen, 1986a, b). 
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C. MESOSCALE EDDIES 



Studies of small-scale (5 to 15 km) eddies observed in conjunction with the ice 
edge front over deep water (Johannessen et al., 1983 and Johannessen et al., 1987) have 
significantly increased the knowledge of the structure of mesoscale eddies. The front 
was observed to have velocities of approximately 10 cm, s (with little vertical shear) and 
a width scale of 10 km. Observations found several mesoscale eddies with length scales 
ranging from 5 to 15 km associated with the front. Like the front, these associated 
eddies also exhibited little vertical shear. Based on considerations from theory and 
laboratory experiments discussed above, Johannessen et al. (1983) reasoned that the 
eddies most likely resulted from barotropic instability of the oceanic front. 
Johannessen et al. (19S7) also established that eddies and meanders are the dominant 
features along the ice edge under moderate wind conditions. Mesoscale eddies are the 
focus of this study. 

D. WIND-FORCED STUDIES 

Numerous studies of wind driven ocean motion in the marginal ice zone have 
been conducted since 1980. Initially these were two-dimensional cross-ice edge 
examinations of upwelling and downwelling produced by along-ice edge winds (Roed 
and O'Brien, 1983; Ilakkinen, 1986a). Because the coupling is stronger under ice, 
Ekman transport is larger under ice. Roed and O'Brien ( 1983) describe the relationship 
between the wind and the ice edge as “ . . . winds which blow along the ice edge and to 
the left when facing the ice which favors upwelling.” Hakkinen (1986b) also showed 
that along ice edge. 10 m s winds produce upwelling (downwelling) when the ice edge is 
to the right (left) of wind direction. 

Recently, Smith et al. (1987) have studied the effect of along-ice edge winds on 
existing ocean eddies. Their results show that for constant 10 m s winds, the ice 
responds largely to the wind and, to a lesser extent, the ocean eddy. They also find, 
however, that upper ocean anticyclones (cyclones) are rapidly destroyed by upwelling 
(downwelling) favorable winds. None of these studies have considered the effects of 
lighter winds and reversals in the wind direction. 

Simulations using constant wind directions are not altogether realistic. Figure 1.2 
and Figure 1.3 show the wind pattern observed by Johannessen et al. (1987) and 
M orison et al. (19S7), respectively, during MIZEX 84. The two observations show 
that there were periods when the wind shifted regularly, and also there were extended 
periods with light winds. 
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Figure 1.2 Time Series of Wind {from Johannessen et al. (1987)). 
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Johannessen et al. (19S7) provide a schematic of the interaction of an open 
ocean eddy with the ice edge (Figure 1.4). This eddy (Johannessen et al. (1987), eddy 
(El 3)) was observed for about eleven days. Figure 1.2 shows that for the first three 
days, steady winds along the ice were observed, followed bv periods of lighter winds, 
varying in direction for periods of about one day. The last days of observation of this 
eddy were dominated by strong winds from the south, forcing the ice away from the 
eddy. 

E. PURPOSE OF THIS STUDY 

The objective of this study is to understand the interaction of mesoscale, and 
largely barotropic, eddies with a marginal ice zone through the use of a nonlinear two- 
layer regional ocean model with temporally varying winds and a coupled ice model. 
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Figure 1.4 Open Ocean Eddy Interacting with Ice Edge. 



From Johannessen et al. (1987). Where AW is Atlantic Water, MAW is MIZ Atlantic 
Water, ASW is Arctic Slope Water, and MASW is MIZ Arctic Slope Water. 
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II. NUMERICAL TECHNIQUE AND MODEL PARAMETERS 



A. OCEAN MODEL 

1. Ocean Model Equations 

Experiments are performed using a numerical model. The model is based on 
primitive equations, has two layers, and uses a semi-implicit numerical scheme. 
Motion in each layer is governed by a momentum equation: 

dv-jd t + (V-V- + V;-V) V; + k X fV; = 

1 1 1^1 1 (eqn 2.1) 

-h-7P- + A h V 2 Vj + (6 il /p 1 )((l-A)T aw + T IW ) 

and a continuity equation: 

ahj/dt + V ' Vj = 0 (eqn 2.2) 

for layer (i= 1 upper and i= 2 lower), thickness hj, transports Vj, and velocities Vj. The 
fluid is hydrostatic, Boussinesq, and the fluid density (P|) in each immiscible layer is 
fixed. Subgrid scale turbulent eddy dissipation processes are represented by a 

horizontal Laplacian operator on transport. As this only approximates a complex 
subgrid scale dissipation process, its coefficient, A^, has been chosen small (10 m 2 /s). 
The fluid is thus relatively inviscid. The notation used in all equations is defined in the 
Appendix. The above scheme has been used in numerous mid-latitude ocean mesoscale 
circulation studies (Hurlburt and Thompson (1980, 1982); Smith and O'Brien (1983)) 
where more thorough presentations show how energy and enstrophy are conserved in 
the absence of dissipation. Smith and Reid (1982) compared test cases to linear 
analytic solutions to verify the Rossby dispersion characteristics of the model. 

2. Boundary and Initial Conditions 

A rectangular (55 x 40 km) finite difference gridded domain is used. There are 
55 grid points in the x-direction and 44 grid points in the y-direction. The domain is 
rotated counter-clockwise approximately 30 degrees relative to a reference latitude of 
80 °N so that the x-axis is aligned with the East Greenland Current, approximately 
parallel to the east coast of Greenland. The initial ocean state consists of a Gaussian 
height field in each layer, which is centered in the basin: 
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(eqn 2.3) 



h . _ Ao e(-(x 2 + y 2 )/2L 2 )) 

L is the e-folding width scale for the eddy ( = 5 km). The depth of thermocline in the 
East Greenland Current ranges from 25 to 100 m. All experiments for this study are 
initialized with the upper layer mean thickness of 50 m and the lower layer mean 
thickness of 4000 m. The resulting first internal Rossby radius of deformation (R^) is 
equal to 5 km, consistent with the observations of Johannessen el al. (1983). 
Nondimensional eddy size, y=L/R ( j, is equal to 1. The amplitude (A Q ) of the 
Gaussian distribution was chosen so that a maximum velocity would be on the order of 
10 cm/s in each layer. Thus, the Rossby number ( v max /fL) for such a flow is 0.15; 
therefore, the nonlinear terms are expected to be important parameters in the 
momentum equations. The ocean velocity field is thus initially in gradient balance and 
constant with depth (Uj = U 2 ). The north and south boundaries of the model domain 
are no-slip walls, where both tangential flow and normal flow are set equal to zero. 
Open radiation is allowed at the east and west boundaries of the model domain to 
allow wind-driven inflow and outflow. Simulations were integrated for short periods, 
typically 4 days. The sense of ocean eddy rotation (cyclonic vs. anticyclonic), wind 
direction, wind speed, and the duration of wind shifts were varied. Although the eddies 
are initialized with barotropic structure, baroclinic structure associated with wind- 
driven layer interface motion is free to evolve. 

B. ICE MODEL 

1. Ice Model Equations 

Motion in the ice is governed by momentum equations: 
dxxjdx + udu/<5x + vdu/dy - fv -(A/m)(t x a * - t x ‘ w ) -g*5(hj + h 2 )/dx ( ec l n 

dvjdx + ud\/dx + \d\jd y = -fu -(A/m)(ty a *- ty* w ) -g^^hj + h 2 )/<?y (eqn 2.5) 

and by continuity equations: 

dAjdx + <)(Au)/<3x + <3(Av) jd y = A a V^A ( ec l n 
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(eqn 2.7) 



dm/dt + d(mu)/dx + d(mv) }dy = A m V 2 m 

for ice concentration (A) and ice mass (m). As in previous ice model studies, a 
Laplacian damping term has been included in the continuity equations for A and m. 
The ice equations contain a pressure force associated with sea surface slope. Most 
mesoscale modeling studies of the MIZ have neglected this term; Smith et al. (1987) 
show that this term is important in the absence of wind forcing. No internal ice stress 
was used, which is suitable for concentrations less than 85%. Thermodynamic effects 
are not considered in this study. The ice is coupled to the ocean through an ice-water 
interfacial stress, T IW , for ice (u) and ocean (u^,) velocity vectors: 

tlW= Pi c iw( u ' u w)l u * u wl (eqn 2.8) 

and similarly the ice is coupled to the air through the air-ice interfacial stress, T ai , for 
air (u a ) and ice (u; ce ) velocity vectors: 

t ai = P a c a i( u a- u ice)l u a- u ic e l ^ n 2 ' 9 ) 

All constant values (p-, p a , cj w , c a j) are chosen following Hakkinen (1986a) and are 
given in the appendix. 

2. Ice Thickness and Concentration 

The ice thickness distribution, D, is initially specified to be 2 meters thick, but 
is then allowed to vary according to 

D = m/(PjA) (eqn 2.10) 

Ice concentration is initially a concentrated band in the middle of the domain. The ice 
band is a linear function of the y-coordinate, varying from 0.05 at the edge of the band 
nearest the open ocean to 0.75 at the edge of the ice band nearest the pack ice. The 
ice is initially at rest. The difference in value between adjacent contour lines of ice 
concentration found in the model output graphics is 0.05. 
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C. PRELIMINARY CONSIDERATIONS AND EXPECTED MODEL 
RESPONSE 

Previous studies of isolated eddies have examined in detail the propagation and 
decay of isolated eddies in quiescent P-plane backgrounds (McWilliams and Flierl, 
1979; Mied and Lindemann, 1979). At mid-latitudes, isolated eddies have significant 
propagation tendencies associated with the planetary vorticity gradient (P) and can 
decay through Rossby wave radiation and frictional processes. These studies showed 
that westward propagation speed is limited by the Rossby long wave speed: 

c max ~ *P ^d^ ^ eC l n 2-11) 

The value of c max at low latitudes is approximately 4 km/day; however, at 80 °N, 
c max * s on or< ^ er -01 km/day. The decay time associated with Rossby wave 
radiation (Smith and Reid, 1982) is substantially longer at high latitudes. Simulations 
reported here focus on the short time interaction (less than 1 week) of isolated eddies 
and ice. It is anticipated that eddy evolution and decay for these short time cases are 
related wind or to ice processes and not to the aforementioned processes, which occur 
on longer time scales. 

The previous wind forced marginal ice zone ocean models give an indication of 
the expected Ekman response in the ocean and ice edge. Hakkinen (1986a, b) shows 
that 10 m/s along-ice-edge winds give an upwelled or downwelled interface amplitude 
of ~10 m and an ice edge jet velocity of 9 cm/s after a 10-15 day spin-up period. 
Additionally, ice banding simulations require at least one wind reversal to occur. As 
simulations here are integrated for shorter periods with weaker winds, weak ice edge jet 
velocities and no ice banding are expected. 
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III. EXPERIMENTS 



A. WIND FORCING DIRECTIONS 

An understanding of the effect of wind direction on ice motion is the goal of the 
following experiments. Variable in the experiments are the wind magnitude, direction, 
and duration for a given direction. Values for wind direction are measured 
counterclockwise in degrees from the positive x-axis of the model domain. Since the 
model domain is rotated 30° counterclockwise, a model input wind direction of 025° 
means that on the globe, the wind would be generally from the southwest. This 
correlates to a downwelling favorable wind condition when viewed relative to the initial 
concentrated ice band. Table 1 lists the wind directions used in simulations. 





TABLE 1 




WIND DIRECTIONS USED IN SIMULATIONS 


RELATIVE 


MODEL 


TRUE 


WIND 


INPUT WIND 


WIND 


Upwelling Favorable 


205° 


235° 


Downwelling Favorable 


025° 


055° 


On-ice 


115° 


145° 


Olf-ice 


295° 


325° 



A model input wind direction of 205° would be a true wind generally from the 
northeast; this correlates to an upwelling favorable wind. Winds perpendicular to the 
downwelling or upwelling wind directions are considered to be generally on-ice and off- 
ice winds. Along-ice-edge wind forcing (0° and 180°) induces a cross-ice-edge Ekman 
drift in the ice. The 025° and 205° angles chosen by Smith et al. (1987) are intended 
to minimize (but not eliminate) the cross-ice-edge motion, thereby maintaining the 
concentrated ice band near mid-domain and thus reducing boundary effects. The most 
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important aspect of these experiments is not the true direction of the wind, but the 
relative direction of the wind in relation to the concentrated ice band. 

B. PRELIMINARY EXPERIMENTS 

Before proceeding to the more realistic wind-driven simulations, the response of 
the ice edge to ocean eddies, without wind forcing, is considered first. 

I. Simulations with No Wind 

a. Cyclonic Eddy Initialization 

Figure 3.1 shows the evolution of the concentrated ice band under the 
influence of a cyclonic ocean eddy without wind forcing. By day one, the eddy is 
affecting the ice band in such a way that the ice band contains a sinuous pattern that 
by day two is more pronounced. By day three the ice clearly exhibits a circular, 
cyclonicallv turning eddy, which is completely distinct from the original ice band by 
day four. 

b. Anticyclonic eddy initialization 

As would be intuitively expected, the anticyclonic case (Figure 3.2) with no 
wind forcing exhibits nearly the same result as above, except that the sense of rotation 
is opposite from the cyclonic case. 

The results indicate that the ice velocity evolves rapidly to the ocean 
velocity and is simply advected radially. These results are consistent with Smith et al. 
(1987), where a broader ice concentration band was considered. Smith et al. (1987) 
show that the ice momentum equation in no wind simulations equilibrates to a 
gradient balance. For the cyclone case, a radially inward pressure force in the ice 
associated with the ocean surface slope is balanced by the radially outward Coriolis 
and centrifugal forces. Likewise, for an anticyclone case, the radially outward pressure 
gradient and centrifugal forces are balanced by the imvard Coriolis force. The ice 
momentum balance was shown to equilibrate rapidly in less than one day. 

C. BASE CASE 

Numerous experiments were conducted using the model previously discussed. To 
assess the model output variations resulting from parameter changes, one simulation is 
chosen to be a base case. This simulation was selected as the reference case because of 
its resemblance to the observations of Johannessen et al. (1987). A summary of the 
model parameters chosen for the base case is listed in Table 2. 
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Figure 3.1 Preliminary Experiment (Cyclone). 



Evolution of (a) ice concentration, (b) ocean upper laver potential vorticitv, (c) ocean 
lower layer potential vorticitv, (d) ocean surlace height anomalv, ancf (e) ocean 
interface height anomalv. Contour interval for ice concentration is 0.05. 
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Figure 3.2 Preliminary Experiment (Anticyclone). 



Evolution of (a) ice . concentration, (b) ocean upper layer potential vorticitv. (c) ocean 
lower layer potential vorticity, (d) ocean surface height anomalv. and' (e) ocean 
interface height anomaly. Contour interval for ice concentration is 0.1)5. 
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